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Abstract 
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1 Introduction 



In credit derivatives trading, one important question is how the market compensates investors for 
bearing credit risk. A number of studies [H [5l [TSl [27] have examined analytically and empirically 
the structure of default risk premia inferred from the market prices of corporate bonds, credit 
default swaps, and multi-name credit derivatives. A major risk premium component is the mark- 
to-market risk premium which accounts for the fluctuations in default risk. Under reduced-form 
models of credit risk [181 [2H1 EI] , this is connected with a drift change of the state variable diffusion 
driving the default intensity. In addition, there is the event risk premium (or jump-to-default risk 
premium) that compensates for the uncertain timing of the default event, and is measured by the 
ratio of the risk- neutral intensity to the historical intensity (see [H [57] ) . 

From standard no-arbitrage pricing theory, risk premia specification is inherently tied to the 
selection of risk-neutral pricing measures. A typical buy-side investor (e.g. hedge fund manager 
or proprietary trader) would identify trading opportunities by looking for mispriced contracts in 
the market. This can be interpreted as selecting a pricing measure to reflect her view on credit 
risk evolution and the required risk premia. As a result, the investor's pricing measure may differ 
from that represented by the prevailing market prices. In a related study, Leung and Ludkovski 
[33] showed that such a price discrepancy would also arise from pricing under marginal utility. 

Price discrepancy is also important for investors with credit-sensitive positions who may need 
to control risk exposure through liquidation. The central issue lies in the timing of liquidation 
as investors have the option to sell at the current market price or wait for a later opportunity. 
The optimal strategy, as we will study, depends on the sources of risks, risk premia, as well as 
derivative payoffs. 

This paper tackles the optimal liquidation problem on two fronts. First, we provide a gen- 
eral mathematical framework for price discrepancy between the market and investors under an 
intensity-based credit risk model. Second, we derive and analyze the optimal stopping problem 
corresponding to the liquidation of credit derivatives under price discrepancy. 

In order to measure the benefit of optimally timing to sell as opposed to immediate liquidation, 
we employ the concept of delayed liquidation premium. It turns out to be a very useful tool for 
analyzing the optimal stopping problem. The intuition is that the investor should wait as long as 
the delayed liquidation premium is strictly positive. Applying martingale arguments, we deduce 
the scenarios where immediate or delayed liquidation is optimal (see Theorem 3.4). Moreover, 
through its probabilistic representation, the delayed liquidation premium reveals the roles of risk 
premia in the liquidation timing. Under a Markovian credit risk model, the optimal timing 
is characterized by a liquidation boundary solved from a variational inequality. For numerical 
illustration, we provide a series of examples where the default intensity and interest rate follow 
Cox-Ingersoll-Ross (CIR) or Ornstein-Uhlenbeck (OU) processes. 

Our study also considers the connection between different risk-neutral pricing measures (or 
equivalent martingale measures) in incomplete markets. Well-known examples of candidate pric- 
ing measures that are consistent with the no-arbitrage principle include the minimal martingale 
measure [21], the minimal entropy martingale measure |22l 123] . and the g-optimal martingale 
measure \25\ [26] . The investor's selection of various pricing measures may also be interpreted via 
marginal utility indifference valuation (see, among others, jlll 1331 [M]). 

For many parametric credit risk models, the market pricing measures and risk premia can 
be calibrated given sufficient market data of credit derivatives. For instance, Berndt et al. [5] 
estimated default risk premia from credit default swap (CDS) rates and Moodys KMV expected 
default frequency (EDF) measure. For CDO tranche spreads, Cont and Minca [9j constructed a 
pricing measure and default intensity based on entropy minimization. In this paper, we focus on 
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investigating the impact of pricing measure on the investor's hquidation timing for various credit 
derivatives, including defaultable bonds, CDSs, as well as, multi-name credit derivatives. 

In recent literature, a number of models have been proposed to incorporate the idea of mispric- 
ing into optimal investment. Cornell et al. [TO] studied portfolio optimization based on perceived 
mispricing from the investor's strong belief in the stock price distribution. Ekstrom et al. |19] 
investigated the optimal liquidation of a call spread when the investor's belief on the volatility 
differs from the implied volatility. On the other hand, the problem of optimal stock liquidation 
involving price impacts has been studied in [H 1391 HOj . among others. 

Our work is closest in spirit to [32j where the delayed purchase premium concept was used to 
analyze the optimal timing to purchase equity European and American options under a stochastic 
volatility model and a defaultable stock model. In contrast, the current paper addresses the 
optimal timing to liquidate various credit derivatives. In particular, we adopt a multi-factor 
intensity-based default risk model for single-name credit derivatives, and a self-exciting top-down 
model for a credit default index swap. As a natural extension, we also investigate the optimal 
timing to buy and sell a credit derivative, with or without short-sale constraint, and provide 
numerical illustration of the the optimal buy-and-sell strategy. 

The rest of the paper is organized as follows. In Section 2, we present the mathematical model 
for price discrepancy and formulate the optimal liquidation problem under a general intensity- 
based credit risk model. In Section 3, we study the problem within a Markovian market and 
characterize the optimal liquidation strategy for a general defaultable claim. In Section 4, we 
apply our analysis to a number of single-name credit derivatives, e.g. defaultable bonds and 
credit default swaps (CDS). In Section 5, we discuss the optimal liquidation of credit default 
index swap. In Section 6, we examine the optimal buy-and-sell strategy for defaultable claims. 
Section 7 concludes the paper and suggests directions for future research. 

2 Problem Formulation 

This section provides the mathematical formulation of price discrepancy and the optimal liq- 
uidation of credit derivatives under an intensity-based credit risk model. We fix a probability 
space {Cl,G,F), where F is the historical measure, and denote T as the maturity of derivatives 
in question. There is a stochastic risk-free interest rate process {rt)o<t<T- The default arrival is 
described by the first jump of a doubly-stochastic Poisson process. Precisely, assuming a default 
intensity process {\t)o<t<T-, we define the default time by 



generated by N (see e.g. ^li Chap. 5]). 
2.1 Price Discrepancy 

By standard no-arbitrage pricing theory, the market price of a defaultable claim, denoted by 
{Pt)o<t<T, is computed from a conditional expectation of discounted payoff under the market 
risk- neutral (or equivalent martingale) pricing measure Q ~ P. In many parametric credit risk 
models, the market pricing measure Q is related to the historical measure P via the default 





(2.1) 
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risk premia (see Section 3.1 below). We assume the standard hypothesis (H) that every F- local 
martingale is a G- local martingale holds under Q (see [7, §8.3]). 

We can describe a general defaultable claim by the quadruple {Y, A, R,Td), where Y £ Ft is 
the terminal payoff if the defaultable claim survives at T, (^t)o<t<T is a F-adapted continuous 
process of finite variation with Aq = Q representing the promised dividends until maturity or 
default, and {Rt)o<t<T is a F-predictable process representing the recovery payoff paid at default. 
Similar notations are used by Bielecki et al. where the following integrability conditions are 
assumed: 

E^{|e-io'""'^^y|} < oo, E^{\ [ e- ^0 '^'^^{1 - Nu)dAJ} <oo, and 

J{Q,T] 

E^{|e-i'o''^^'^''*i?,,AT|} < oo. (2.2) 
For a defaultable claim {Y, A, R,T(i), the associated cash flow process {Dt)o<t<T is defined by 

A := l"l{r,>T}l{t>T} + / {l-Nu)dAu+ j RudNu. (2.3) 

J(0,t^T] J{0,tAT] 

Then, the (cumulative) market price process {Pt)o<t<T is given by the conditional expectation 
under the market pricing measure Q : 

J(Q,T] 

One simple example is the zero-coupon zero-recovery defaultable bond (l,0,0,rrf), whose market 

price is simply P* = E^{e- ''-'^''l{r^>T}\Gt] ■ 

When a perfect replication is unavailable, the market is incomplete and there exist different 
risk-neutral pricing measures that give different no-arbitrage prices for the same defaultable claim. 
Mathematically, this amounts to assigning a different risk-neutral pricing measure Q ~ Q. The 
investor's reference price process {Pt)o<t<T is given by the conditional expectation under investor's 
risk-neutral pricing measure Q : 

Pt:=E'i{[ e-Ii"'^''^dD^\gt}, (2.5) 

i(0,T] 

whose discounted price process (e~ ■^o''''^''Pj)o<t<T is a (Q, G)-martingale. We assume that the 
standard hypothesis (H) also holds under Q. 



2.2 Optimal Stopping & Delayed Liquidation Premium 

A defaultable claim holder can sell her position at the prevailing market price. If she completely 
agrees with the market price, then she will be indifferent to sell at any time. Under price discrep- 
ancy, however, there is a timing option embedded in the optimal liquidation problem. Precisely, 
in order to maximize the expected spread between the investor's price and the market price, the 
holder solves the optimal stopping problem: 

Jt := ess sup E^{e" "''^'^"{Pr - Pr)\gt}, < t < T, (2.6) 
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where Tt^r is the set of G-stopping times taking values in [t,T]. Using repeated conditioning, we 



decompose (2.Q) to Jt = Vt — Pt, where 



Vt := esssup E^{e- -^^ '^^'^'^ PrlGt} ■ (2.7) 



Hence, maximizing the price spread in (2.6) is equivalent to maximizing the expected discounted 



future market value Pt under the investor's measure Q in (2.7). 

The selection of the risk-neutral pricing measure Q can be based on the investor's hedging 
criterion or risk preferences. For instance, dynamic hedging under a quadratic criterion amounts to 
pricing under the well-known minimal martingale measure developed by Follmer and Schweizer 
|21j . On the other hand, different risk-neutral pricing measures may also arise from marginal 
utility indifference pricing. In the cases of exponential and power utilities, this pricing mechanism 
will lead the investor to select the minimal entropy martingale measure (MEMM) (see [33]) and 
the (/-optimal martingale measure (see [25]). 

Lemma 2.1. ForO<t<T, we have Vt> PtV Pt- Also, Vr^ = Pr^ = Pr^ at default. 



Proof. Since t = t and t = T are candidate liquidation times, we conclude from (2.7) that 
Vt>PtV Pt. Also, we observe from that Pt = /^g e" -A" '''''^"dDu = Pt for t>TdAT. This 
implies that 

= esssup E^{e~^-d'^^'^''Pr\gr^} = esssup E^{e" ^-<i ''"'^''P^ig^,} = = Pr^. (2.8) 



□ 



The last equation means that price discrepancy vanishes when the default event is observed or 
when the contract expires. This is also realistic since the market will no longer be liquid afterward. 
If the defaultable claim is underpriced by the market at all times, that is, Pt < -Pt, Vi < T, 



then we infer from (2.6) that Jt = 0. This can be achieved at r* = T since price discrepancy must 



vanish at maturity, i.e. Pt = Pt- In turn, this implies that 

Vt = E^{e--ft^''-'^''PT\gt} = E^{e'-fi'''-'^''PT\gt} = Pt- 
In this case, there is no benefit to liquidate before maturity T. 



According to (2.7), the optimal liquidation timing directly depends on the investor's pricing 
measure Q as well as the market pricing measure Q (via the market price P). Specifically, we 
observe that the discounted market price {e~ ■l^o^"'^'" P^^Q^f^rp is a (Q, G)-martingale, but generally 
not a (Q, G)-martingale. If the discounted market price is a (Q, G)-supermartingale, then it is 
optimal to sell the claim immediately. If the discounted market price turns out to be a (Q, G)- 
submartingale, then it is optimal to delay the liquidation until maturity T. Besides these two 
scenarios, the optimal liquidation strategy may be non-trivial. 

To quantify the value of optimally waiting to sell, we define the delayed liquidation premium: 

Lt := Vt-Pt> 0. (2.9) 

It is often more intuitive to study the optimal liquidation timing in terms of the premium L. 
Indeed, standard optimal stopping theory [291 Appendix D] suggests that the optimal stopping 
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time r* for (2.7) is the first time the process V reaches the reward P, namely, 

T* = inf{t <u<T :Vu = Pu} = inf{t < n < T : L„ = 0}. (2.10) 



The last equation, which follows directly from definition (2.9), implies that the investor will 



liquidate as soon as the delayed liquidation premium vanishes. Moreover, we observe from (2.8) 



and (2.10) that t* < r^. 



3 Optimal Liquidation under Markovian Credit Risk Models 

We proceed to analyze the optimal liquidation problem under a general class of Markovian credit 
risk models. The description of various pricing measures will involve the mark-to-market risk pre- 
mium and event risk premium, which are crucial in the characterization of the optimal liquidation 



strategy (see Theorem 3.4). 



3.1 Pricing Measures and Default Risk Premia 

We consider a n-dimensional Markovian state vector process X that drives the interest rate 
rt = r(t, X() and default intensity A* = A(t, Xt) for some positive measurable functions r(-, •) and 
A(-, •). Denote by F the filtration generated by X. We also assume a Markovian payoff structure 
for the defaultable claim {Y,A,R,Td) with Y = y(Xj'), At = Jq q{u,X.u)du, and Rt = R{t,X.t) 



for some measurable functions Y(-), q{-, •), and R{-, •) satisfying integrability conditions (2.2) 
Under the historical measure P, the state vector process X satisfies the SDE 

dXt = a{t,Xt)dt + ^{t,Xt)dWf, (3.1) 

where W"" is a m-dimensional P-Brownian motion, a is the deterministic drift function, and S is 
the n by m deterministic volatility function. Standard Lipschitz and growth conditions [301 §5.2] 



are assumed to guarantee a unique solution to (3.1). 



Next, we consider the market pricing measure Q P. To this end, we define the Radon- 
Nikodym density process (Z^''^)o<t<T by 

= ^1^* = ^( - 0^'"-W^),£:((/z - l)M^)^ , (3.2) 
where the Doleans-Dade exponentials are defined by 

£{ - 0^'^. W^), := exp( - i j^lcl^rfdu - ^ V^-dW^) , (3.3) 

£{{fi - 1)M^)^ := exp(^j\og{fiu-)dNu - j\l - N^Kfiu - l)A„du^, (3.4) 

and Mf := Nt — Jq(1 — Ny)\udu is the compensated (P, G)-martingale associated with N. Here, 
{4>t'^)o<t<T and {fit)o<t<T are adapted processes satisfying \\cf)2'^pdu < oo, ^ > 0, and 
Jq fJ'u^udu < oo (see Theorem 4.8 of [41j). 

The process (f)'^'^ is commonly referred to as the mark-to-market risk premium (see [4j), which 
is assumed herein to be Markovian of the form cf)^'^(t, X^). The process fx is referred to as event 
risk premium (see [11[27]), which captures the compensation from the uncertain timing of default. 
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The Q-default intensity, denoted by A, is related to P-intensity via \t = fJ-tXt- Here, we also 
assume /i to be Markovian of the form /i(t,Xt) = A(t, Xj)/A(t, Xt). 

By multi-dimensional Girsanov Theorem, it follows that wj^ := Wf + cpl^'^du is a m- 
dimensional Q-Brownian motion, and M^^ := Nt — /o(l — Nu)fiu^udu is a (Q, G)-martingale. 
Consequently, the Q-dynamics of X are given by 

dXt = b{t, Xt)dt + S(t, Xt)(iWp, (3.5) 

where 6(t,Xt) := a{t,Xt) - S(t, Xt)</.^'P(t, X*). 

Similarly, the investor's pricing measure Q is related to the historical measure P through the 
investor's Markovian risk premium functions cf)^'^{t,x) and p,{t,x). Precisely, the measure Q is 
defined by the density process Z'j^''^ = £(^ — (f/^'^ -W^) ^ £ (^{fi — 1)M^^^. By a change of measure, 

the drift of X under Q is modified to 6(t,Xt) := a{t,Xt) - T,{t,Xt)^^^^{t,Xt). 
Then, the EMMs Q and Q are related by the Radon-Nikodym derivative: 

^Q,Q ^ ^1^^ ^ ^( _ 0Q,Q.wQ)^f ((^ - 1)M^)^ , (3.6) 

where the Doleans-Dade exponentials are defined by 

£(-cl^^'^.W^)-.= exp(-U\\ct>2^^fdu- /V?'^-dW?Y (3.7) 



Jo 

t r, ft 



£{{^-l)M^) ■.= exJ [ log(^)diV„- / (l-N^){^-l)X^du). (3.8) 

We observe that (^^'"^ = cfy^'^ — cfy^'^ from the decomposition: 

cfyf'^dt = dwf - dWf = (dwf - dWf ) - {dWf - dWj) = (0p - cf>f'^)dt. (3.9) 

Therefore, we can interpret 0*^'"^ as the incremental mark-to-market risk premium assigned by 
the investor relative to the market. On the other hand, the discrepancy in event risk premia is 
accounted for in the second Doleans-Dade exponential (3.8). 

Example 3.1. The OU Model. Suppose (r, A) = X, following the OU dynamics: 

drA _ ( kriOr -rt)\,^^( Or A / dW]^\ ^ 

d\t) - Ua(^a - A,); + Uap -aTt^j ydw-"/) ' ^^-'^^ 

with constant parameters kj., Or, k\, Ox > 0. Here, kr, kx parameterize the speed of mean reversion, 
and Or, Ox represent the long-term means (see |41l §7.1.1]). Assuming a constant event risk 
premium fi by the market, the Q-intensity is specified by Xt = iJ.Xt and the pair (r. A) satisfies 
SDEs: 

dn\ ^ ( - n)\ ^ (ar OA / dwi^^\ 

with constants Hr,Or, i^x,Ox > 0. Under the investor's measure Q, the SDEs for rt and Xt = fiXt 
are of the same form with parameters kr,0r,kx. Ox and /2, and is replaced by W^. 



8 



Direct computation yields the relative mark-to-market risk premium: 

CT,. 

1 Kx{ex-~Xt)-iix(ex-Xt) _ p K.r(er-rt)-K.r{er-n) 
The upper term is the incremental risk premium for the interest rate while the bottom term 



reflects the discrepancy in the default risk premia (see (3.9)). 

Example 3.2. The CIR Model. Let X = follow the multifactor CIR model [HI 

§7.2]: 



dXl = ki{9i - Xl)dt + ai^Xl dW'/, (3.12) 

where W^"^ are mutually independent P-Brownian motions and Kj, 6i, ai > 0, i = l,...,n 
satisfy Feller condition 2ki9i > af. The interest rate r and historical default intensity A are 



non- negative linear combinations of X* with constant weights wl,w^ > 0, namely, rt = w^Xl 

n 

and At = ''^^wfxl. Under measure Q, X^ satisfies the SDE: 
1=1 

dXl = Ki{e^- Xi )dt + ai dW'l^ , (3.13) 

with new mean reversion speed Kj and long-run mean Qi. 

Under the investor's measure Q, the SDE for the state vector is of the same form with new 
parameters Kj, Qi. The associated relative mark-to-market risk premium has following structure: 

,Q,Q) _ i^iiQj — XI) — ki{6i — XI) 



The event risk premia {/i, jl) are assigned via A^ = fiXt under Q and Aj = flXt under Q respectively. 
Remark 3.3. The current framework can be readily generalized to the situation where the investor 
needs to assume an alternative historical measure P. The resulting risk premium 0*^'"^ will have 
a third decomposition component (f/'''^ , reflecting the difference in historical dynamics. 

For any defaultable claim (Y, A, R,T(i), the ex-dividend pre-default market price is given by 
(see [3 Prop. 8.2.1]) 

C{t,Xt) = E«2{e-/*^(^"+^'')^^y(XT) + j\-^^^'-^+^^^''^{\uR{u,y.u) + q{u,^u))du\Ft}. (3.14) 
The associated cumulative price is related to the pre-default price via 

Pt = {l- Nt)C{t,Xt) + [\l - Nu)q{u,Xu)ef>^''^du + f R{u,X^)ef>^''^dN^. (3.15) 

Jo J{0,t] 

The price function C(t,x) can be determined by solving the PDF: 

{dC 
— (t,x)+£b,AC(t,x)+A(t,x)ii(t,x)+g(t,x) =0, (t,x) G [0,r) xM", ^^^^^^ 
c(r,x) = y(x), xeM", 
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where Lx is the operator defined by 

Cb,xf = + 2 (^(^'''^^^^'''^^^^^ S^ - (^(*''') + ^(*'^))/- (3.17) 

i=l * *ii = l * 

The computation is similar for the investor's price under Q. 

3.2 Delayed Liquidation Premium and Optimal Timing 



Next, we analyze the optimal liquidation problem V defined in (2.7) for the general defaultable 
claim under the current Markovian setting. 

Theorem 3.4. For a general defaultable claim {Y, A, R, r^) under the Markovian credit risk model, 
the delayed liquidation premium admits the probabilistic representation: 

Lt = l{t<r,}esssnp re-/"(^'"+^")'^^G(n,X„)dn|J-t}, (3.18) 

where G : [0,T] x W ^ R is defined by 
G(t,x) = -(V,C(t,x))^S(i,x)(/)«'^(t,x) + (i?(t,x) - C(i,x))(/i(t,x) -/i(t,x))A(t,x). (3.19) 

If G{t,x.) > V(t,x), then it is optimal to delay the liquidation till maturity T. 
If G{t,x.) < V(t,x), then it is optimal to sell immediately. 

Proof. First, we look at the Q-dynamics of discounted market price (e~ Jt Pu)t<u<T- Applying 
Corollary 2.2 of [6j, for t < u < T, 

die-ft^'r^^^Pu) = e-/"^'"^l(i2, - Cu)dM2 + (1 - 7V„)(V,C„)^S„dWa (3.20) 
= e-i'^"'^'^((l - Nu)G{u,Xu)du+{l - Nu){VxGuf^udW2 + {Ru - C„) 



where G is defined in (3.19), and M'^ is the compensated (Q, G)-martingale for N. Consequently, 



Lt = ess sup Hi - Nu)e~f>-'^''G{u,J(.u)du\gt}, 



where (3.18) follows from the change of filtration technique [71 §5.1.1]. If G > 0, then the integrand 



in (3.18) is positive a.s. and therefore the largest possible stopping time T is optimal. If G < 0, 



then r* = t is optimal and Lt = a.s. □ 

The drift function G has two components explicitly depending on cf)^'^ and fl — fj,. If 
(f)'^'^(t,x) = V(t,x), that is, the investor and market agree on the mark-to-market risk pre- 
mium, then the sign of G is solely determined by the difference fl — n, since recovery R in general 
is less than the pre-default price C. On the other hand, if /i(t, x) = fj,{t, x) V(t, x), then the second 
term of G vanishes but G still depends on /x through V^C in the first term. 



Theorem 3.4 allows us to conclude the optimal liquidation timing when the drift function is 
of constant sign. In other cases, the optimal liquidation policy may be non-trivial and needs 
to be numerically determined. For this purpose, we write Lt = l|f<^^}L(t, Xt), where L is the 
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(Markovian) pre-default delayed liquidation premium defined by 



L(t,Xt) = esssup / e-i'"(''"+^")'^''G(u, X„)du| J^}. 

T&Tt,T Jt 



(3.21) 



We determine L from the variational inequality: 
dL 



mm 



dt 



(t,x)-£5_3^L(t,x)-G(t,x), L(t,x))=0, (t,x)G[0,r)x 



(3.22) 



where £r t is defined in (3.17), and the terminal condition is L{T,:x.) = 0, for x G 



The investor's optimal timing is characterized by the sell region S and delay region P, namely, 



5 = {(t,x) G [o,r] X 
P = {(t,x) G [0,T] X 



L(t,x) =0}, 
L(f,x) > 0}. 



(3.23) 
(3.24) 



Also, define f* = inf{t < u < T : = 0} . On {f* > t^}, liquidation occurs at since L^^ = 0. 
On {f* < Td}, T* is optimal since when u < f*, = l|„<^^}.Lu > and Lf* = 0. Incorporating 
the observation of r^, the optimal stopping time is r* = f* A r^. 

Hence, given no default by time t and = x, it is optimal to wait at the current time t if 
-L(t, x) > in view of the delay region D in (3.24). This is also intuitive as there is a strictly 



positive premium for delaying liquidation. On the other hand, the sell region S must lie within 
the set G- := {(t,x) : G{t,x) < 0}. To see this, we infer from (3.21) that, for any given point 
(t,x) such that L(t,x) = 0, we must have G{t,x.) < 0. In turn, the delay region D must contain 
the set G+ := {{t, x) : G{t, x) > 0}. From these observations, one can obtain some insights about 
the sell and delay regions by inspecting G(t, x), which is much easier to compute than L{t,x). 
We shall illustrate this numerically in Figures 1-4. 

Lastly, let us consider a special example where the stochastic factor X is absent from the model. 



With reference to (3.14), we set a constant terminal payoff Y, and deterministic recovery R{t) 



and coupon rate q{t). Suppose the investor and market perceive the same deterministic interest 
rate r(t), but possibly different deterministic default intensities, respectively, X{t) = fl{t)X{t) and 
X{t) = n{t)X{t). In this case, the price function C in (3.19) will depend only on t but not on x, 



and there will be no mark-to-market risk premium. Therefore, the first term of drift function in 



(3.19) will vanish. However, the second term remains due to potential discrepancy in event risk 
premium, i.e. fl{t) ^ fJ-it). As a result, the drift function reduces to 

G{t) = {R{t)-c{t)mt)-mm. 

Furthermore, the absence of the stochastic factor X also trivializes the filtration F, and leads the 
investor to optimize over only constant times. The delayed liquidation premium admits the form: 



Lt 



|L(t), where L{t) is a deterministic function given by 



L{t) 



sup 

t<i<T 



(3.25) 



As in Theorem 3.4 if G is always positive (resp. negative) over [t, T], then the optimal time t* = T 
(resp. t* = t). Otherwise, differentiating the integral in (3.25) implies that the deterministic 
candidate times also include the roots of G{t) = 0. Therefore, we select among the candidate 
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times t,T and the roots of G to see which would yield the largest integral value in (3.25). 



4 Application to Single-Name Credit Derivatives 

We proceed to illustrate our analysis for a number of credit derivatives, with an emphasis on how 
risk premia discrepancy affects the optimal liquidation strategies. 



4.1 Defaultable Bonds with Zero Recovery 

Consider a defaultable zero-coupon zero-recovery bond with face value 1 and maturity T. By a 
change of filtration [71 §5.1.1], the market price of the zero-coupon zero-recovery bond is given by 

:= E'^{e-/^"'^n|,^>r||gi} = l|i<,^}EQ{e-/"(^"+^")'^''|J-<} = l{,^,^yC\t,Xt), (4.1) 
where denotes the market pre-default price that solves ( 3.16[ ). Under the general Markovian 



credit risk model in Section 3.1, we can apply Theorem 3.4 with the quadruple (1,0,0, r^) to 



obtain the corresponding drift function. 

Under the OU dynamics in Section 3.1, the pre-default price function C^{t,r,X) is given 
explicitly by fiT, §7.1.1]: 

C^(t, r. A) = ^AiT-t)-BiT-t)r-DiT-t)X^ (4 2) 

where 

B{s) = , D{s) = , (4.3) 

A{s) = [ [la^B^{z) + pnaraxB{z)D{z) + J^V^Z)2(z) - k,.0^P(z) - pKxexD{z)]dz. 

As a result, the drift function G^{t, r, A) admits a separable form: 

G°(t, r, A) = C^{t, r, A) ( B{T - t){Kr - Kr)r + B{T - t)(K^e^ - Rr 



+ [D{T - t){kx - ^x) -{-- 1)]A + tiD{T - t){Kx9x - kxOx) ] ■ (4.4) 

We can draw several insights on the liquidation timing from this drift function. If the market 
and the investor agree on the speed of mean reversion for interest rate, i.e. Kr = k^, then 
G'^(t, r, A)/C'^(t, r, A) is linear in A. Furthermore, if the slope D{T — t){kx — kx) ~ (^ ~ 1) ^-iid 

intercept B{T — t){KrOr — krOr) + fJ-DiT — t){K,xOx — kxOx) are of the same sign, then the optimal 
liquidation strategy must be trivial in view of Theorem |3.4[ In contrast, if the slope and intercept 
differ in signs, the optimal stopping problem may be nontrivial and the sign of the slope determines 
qualitative properties of optimal stopping rules. For instance, suppose the slope is positive. We 
infer that it is optimal for the holder to wait at high default intensity where the corresponding 
and thus delayed liquidation premium are positive. The converse holds if the slope is negative. 

If the investor disagrees with market only on event risk premium, i.e. n ^ fl, then the drift 
function is reduced to G^{t,r,X) = —C^{t,r,X){^ — 1)A, which is of constant sign. This implies 
trivial strategies. If ^ > fl, then G*^ > and it is optimal to delay the liquidation until maturity. 
On the other hand, if /i < /i, then it is optimal to sell immediately. More general specifications 
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of the event risk premium could depend on the state vector and may lead to nontrivial optimal 
stopping rules. Disagreement on mean level 0\ has a similar effect to that of fx. 

If the investor disagrees with market only on speed of mean reversion, i.e. kx 7^ kx, then 
r, A) = C°(t, r, X)D{T - t) [{kx - n\)X + /x6Ia(ka - ^a)] with D{T - t) > before T, where 
the slope and intercept differ in signs. If ka < k^, the slope Rx — nx is positive and it is optimal to 
sell immediately at a low intensity, and thus, a high bond price. The converse holds for nx > ^A- 

We consider a numerical example where the interest rate is constant and the market default 
intensity A is chosen as the state vector X with OU dynamics. W e emp loy the standard implicit 
PSOR algorithm to solve L{t, A) through its variational inequality (3.22 ) over a uniform finite grid 
with Neumann condition applied on the intensity boundary. The market parameters are T = 1, 
fi = 2, Kx = 0.2, 9x = 0.015, r = 0.03, and a = 0.02, which are based on the estimates in OttS]. 



From formula (4.2), we observe a one-to-one correspondence between the market pre-default 
bond price and its default intensity A for any fixed (t, r), namely, 



A 



-log(C°) + A{T -t)- B{T - t)r 
D{T - t) ■ 



(4.5) 



Substituting (4.5) into (3.23) and (3.24), we can characterize the sell region and delay region in 



terms of the observable pre-default market price 

In the left panel of Figure [T| we assume that the investor agrees with the market on all 
parameters, but has a higher speed of mean reversion kx > kx- In this case, the investor tends 
to sell the bond at a high market price, which is consistent with our previous analysis in terms 
of drift function. If the bond price starts below 0.958 at time 0, the optimal liquidation strategy 
for the investor is to hold and sell the bond as soon as the price hits the optimal boundary. If the 
bond price starts above 0.958 at time 0, the optimal liquidation strategy is to sell immediately. In 
the opposite case where kx < kx (see Figure [l](right)), the optimal liquidation strategy is reversed 
- it is optimal to sell at a lower boundary. In each cases, the sell region must lie within where G 
is non-positive, and the straight line defined by G = can be viewed as a linear approximation 
of the optimal liquidation boundary. 

admits closed-form formula |4H §7.2] 



Under the CIR dynamics in Section 



3.1 



where 



y^.w^)Xldv^y^^ = x} = J[Ai{T-t)e 



Bi{T-t)xi 



1=1 



25 g(S,+«;,)s/2 



As a result, the drift function is given by 



(4.6) 



(4.7) 
(4.8) 



GO(t,x) = {[B,{T - t){ki - K,) - (/i - f,)w^]x^ + Bi{T - 0(^1^0^ - kA))]C\t,^) 



i=l 



which is again linear in terms of C^(t,x). 

To illustrate the optimal liquidation strategy, we consider a numerical example where interest 
rate is constant, X=A, and = ^. The benchmark specifications for the market default intensity 
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Figure 1: Optimal liquidation boundary in terms of market pre-default bond price under OU dynamics. We 
take T = 1, r = 0.03, a = 0.02, ^ = /i = 2, and 6'a = = 0.015. Left panel: When ka = 0.2 < 0.3 = ka, 
the optimal boundary increases from 0.958 to 1 over time. Right panel: When ka = 0.3 > 0.2 = ka, the 
optimal boundary increases from 0.927 to 1 over time. The dashed straight line is defined by G = 0, and 
we have G < in both sell regions. 



A in the CIR dynamics are T = 1, fj, = 2, kx = 0.2, 9x = 0.015, r = 0.03, and a = 0.07 based on 
the estimates from |14l I15j . Like in the OU model, we can again express the sell region and delay 
region in terms of the pre-default market price see Figure [2I 



4.2 Recovery of Treasury and Market Value 

Extending the preceding analysis on defaultable bonds, we incorporate two principle ways of 
modeling recovery: the recovery of treasury or market value. 

By the recovery of treasury, we assume that a recovery of c times the value of the equivalent 
default-free bond is paid upon default. Therefore, the market pre-default bond price function is 



C^^{t, x) = (1 - c)C"{t, x) + c(3{t, x), 



(4.9) 



where /3(t,x) := E^^je -^t '^"°'^|Xt = x} is the equivalent default-free bond price. Then, applying 
Theorem |3.4| with the quadruple (1, 0, c/3, r^), we obtain the corresponding drift function: 



G^^it, x) = - {V^C^^it, x))^S(t, x)(/)^'Q(t, x) + (c - 1) {il{t, x) - ^l{t, x)) X{t, Jc)C%t, x). 

(4.10) 

If c = 0, then C^^(t,x) = CO(t,x) and G^^ in (4.10) reduces to the drift function of the zero- 
recovery bond. If c = 1, then C^-^(t,x) = /3{t,x.) is the market price of a default-free bond, and 
risk premium discrepancy may arise only from the interest rate dynamics. 

Here are two examples where the drift function G^^ in (4.10) can be computed explicitly. 

Exam ple 4.1. Under OU model, C^'^ {t,r, X) is computed according to (4.9) with C^{t,r,X) in 
^ and /3(t,r,A) = e^(^-*)-^(^-*>, where A{s) = J^l^a'^B^iz) - Kr9rB{z)] dz and B{s) is 
defined in (4.3). 

Example 4.2. Under the multi-factor CIR model, C^' ^{t, x ) is found again from (4.9), where 
CO(f,x) is given in (4.6), and /3(f,x) is computed from (4.6) with w** = in (4.7) and (4.8). 

As for the recovery of market value, we assume that at default the recovery is c times the 
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Figure 2: Optimal liquidation boundary in terms of market pre-default bond price under CIR dynamics. 
Wc take T = 1, r 0.03, a = 0.07, ^Ji ^ jl ^ 2, andOx ^ h ^ 0.015. Left panel: When ka = 0.2 < 0.3 = ka, 
the optimal boundary increases from 0.948 to 1 over time. Right panel: When k,\ = 0.3 > 0.2 = ka, the 
optimal boundary increases from 0.935 to 1 over time. 

pre-default value C^^l^ . The market pre-default price is given by 

C^^^^(t,Xt) = E^{e-^*'^(''''+(^-^)^-)°'nJ-t}, < t < T. (4.11) 

The corresponding drift function can be obtained by applying the quadruple (1, 0, cC^^^^, Td) to 
Theorem 13.41 



Example 4.3. Under the OU model in Section 3.1, the price function C^^^^ {t,r,X) is given by 



where B{s) is defined in (4.3), 



fl - cKl - e"''^'') 
D{s) = ^ '-, and 



Jo ^ ^ 



Example 4.4. Under the multi- factor CIR model, C^^^^(t,x) admits the same formula as (4.6) 
but with W'^ replaced by (1 — c)w^ in (4.7) and (|4.8[). 



4.3 Optimal Liquidation of CDS 

In this section we consider optimally liquidating a digital CDS position. The investor is a pro- 
tection buyer who pays a fixed premium to the protection seller from time until default or 
maturity T, whichever comes first. The premium rate p^, called the market spread, is specified 
at contract inception. In return, the protection buyer will receive $1 if default occurs at or before 
T. The liquidation of the CDS position at time t can be achieved by entering a CDS contract as 
a protection seller with the same credit reference and same maturity T at the prevailing market 
spread p™. By definition, the prevailing market spread makes the values of two legs equal at 
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time t, i.e. 



-{ l\-^^^'-^^Wl{u<r,}du\gt} = EQ{e-/*^'^'"^n|,<,^<r}|gt}. (4.12) 



If the liquidation occurs at time t, she receives the premium at rate p"^ and pays the premium 
at rate until default or maturity T. If default occurs, then the default payments from both 
CDS contracts will cancel. Considering the resulting expected cash flows and (4.12), the mark- 
to-market value of the CDS is given by 

E^^{ fe-^>^''-{pT-p'^)l{u<r,}du\Gt] =l{t<.,}IE«{ /V/>"+^")'^-(A„-pS^)dtx|J-t} 
Jt Jt 

= :l|i<.,}C^^^(t,Xi). (4.13) 
For CDS, we apply the quadruple (0, — p™, 1, t^) to Theorem |3.4| and obtain the drift function: 

G^^^(t,x) = -(V,C^^^(t,x))^S(t,x)0^'^^(t,x) + (1 - C^^^(t,x))(/i(f,x) - /x(t,x))A(t,x). 

(4.14) 

If there is no discrepancy over mark-to-market risk premium, i.e. cf)^''^{t,x) = 0, then the sign 
of G'-"^'^ is determined by /i(t, x) — /x(t, x) since 

(jCDS < Yrom this we infer that higher event 
risk premium (relative to market) implies delayed liquidation. 

In general, the market pre-default value C^^^ can be solved by PDE (3.16). If the state vector 
X admits OU or CIR dynamics, C*^^*^, and thus G^^^ , is given in closed form, as illustrated in 
the following two examples. 



Example 4.5. Under the OU dynamics, the pre-default value of CDS (see (4.13)) is given by the 
following integral: 



C^^^(t,r,A) 



+ 



'^^^'^-'^gis,u)ds-p^ 



du, 



where G^{t, r, A; u) is given by (4.2) with T = u and 



gis,u) := HKx9x - pf^drcrx 



I _ g-Kr(M-s) 



I _ g-KA(«-s) 



Example 4.6. Under the multi- factor CIR dynamics, the pre-default value of CDS is given by the 
following integral: 



C7^^^(t,x)= / C7°(t,x;u) Y,Lw^{KieiBiiu-t)+Bl{u-t)xi)'] 



Po 



du, 



where C''(t, x;m) is given in (4.6) with T = u and Bi{s) in (4.8). See Chapter 7 of [IT] . 

Example 4.7. For a forward CDS with start date Tq < T, the protection buyer pays premium 
at rate pa from until or maturity T, and receives 1 if G [Ta,r]. By direct computation, 
the pre-default market value is C^^^{t, x; T) — C'^^^{t, x; T^), t < Ta- Consequently, closed-form 



formulas for the drift function are available under OU or CIR dynamics by Examples 4.5 and 4.6 



We consider a numerical example where interest rate is constant and state vector X=A follows 
the CIR dynamics. We assume that the investor agrees with the market on all parameters except 
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the speed of mean reversion for default intensity. In the left panel of Figure [3] with k\ = 0.2 < 
0.3 = k\, the optimal liquidation strategy is to sell as soon as the market CDS value reaches an 
upper boundary. In the case with kx = 0.3 > 0.2 = kx (see Figure [3] (right)), the sell region is 
below the continuation region. 
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Figure 3: Optimal liquidation boundary in terms of market pre-default CDS value under CIR dynamics. 
We take T = I, r ^ 0.03, a = 0.07, = 0.02, fi = fi ^ 2, and 9x = 9x ^ 0.015. Left panel: When 
K\ = 0.2 < 0.3 = k\, liquidation occurs at an upper boundary that decreases from 0.0172 to over 
t e [0, 1]. Right panel: When k\ = 0.3 > 0.2 = k\, the CDS is liquidated at a lower liquidation boundary, 
which decreases from 0.00338 to over time. In both cases, the dashed line defined by G = lies within 
the continuation region. 



Remark 4.8. As a straightforward generalization under our framework, one can replace the unit 
payment at default by i?'^^'^(rrf, X^-^). Then, we can apply the quadruple {0, —p"^, R^^^ ,Td) 
to Theorem 3.4, and obtain the same drift function G^^^ in (4.14) except with 1 replaced by 
ii^^^(t,x). 



4.4 Jump-Diffusion Default Intensity 

We can extend our analysis to incorporate jumps to the stochastic state vector. To illustrate this, 
suppose the default intensity and interest rate are driven by a n-dimensional state vector X' with 
the affine jump-diffusion dynamics: 



dX; = a(t, X'Jdt + S(t, X;)dWf + dJt, 



(4.15) 



where J = ( J^, . . . , J")-^ is a vector of n independent pure jump processes taking values in 
M". Under historical measure P, we assume Markovian jump intensity of the form A(i,X^) = 
(A^(t, X^), . . . , A"(i, XJ))-^ for J. All random jump sizes (^Oij of J ^^^^ independent, and for each 
J*, the associated jump sizes Yl,Y2, ■ ■ ■ have a common probability density function /*. 

The default intensity of defaultable security is given by A(i, X^) for some positive measurable 
function A(-,-), and the default counting process associated with default time r^^ is denoted by 



Nt 



L{t>rd}' 



We denote {Gt)o<t<T to be the full filtration generated by W , J, and r^. 



We define a market pricing measure Q in terms of the mark-to-market risk premium <^^' and 
event risk premium fi, which are Markovian and satisfy f^^'^'pfin < oo and ^yX^du < oo. 
Due to the presence of J, the market measure Q can scale the jump intensity of J by the positive 
Markovian factors SI = S^{t,X.[), with j'^ 5l^K.\^du < oo for i = 1, . . . , n. Also, Q can transform 
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the jump size distribution of J by a function (/i*)i=i,...,n > satisfying h^{y)p{y)dy = 1 for 
i = 1, . . . ,n. 

The Radon-Nikodym derivative is given by 

= £{- £:((;. - l)M^)^Kr, 

where first two Doleans-Dade exponentials are defined in (3.3) and (3.4) respectively, Mf := 
Nt — /q(1 — Nu)Xudu is the compensated P-martingale associated with N, and the last term 



A',"-' 



n 

i=l 



exp 



(0 



[l-6\u,K)h\y))A'{u,K)f{y)dydu) J] {S'(T^,X!^,)h'(Yj)) 



where Tj is the jth jump time of r and A^^^'^ := ^ • 

>i l{T»<t} is the counting process associated 
with J*. 

By Girsanov Theorem, Wp := Wf + (p'^'^du is a Q-Brownian motion, the jump intensity 
of r under Q is Ani,X'J := (5^(t, X'JA^(t, X'J, and the jump size pdf of J* under Q is f{y) := 
h^{y)p{y). The Q-dynamics of state vector X' is given by 



dX't = b{t, X.[)dt + X.[)dWf + dJt. 



(4.16) 



Also incorporating the event risk premium, the Q-default intensity is A(t, X^) := Xj)A(t, Xj). 

Under the investor's measure Q, we replace b with b, cj)'^''^ with cf)'^'^ , and W*^ with W'^ for 
the dynamics of X'. For each J\ the Q-intensity is denoted by A*(t,X'j) := Xj)A*(t, XJ), 
and the jump size pdf under Q is := h^{y)P{y). With investor's event risk premium /i, the 

default intensity under Q is A(f,X'J := Xj)A(t, Xj). 

The two pricing measures Q and Q are related by the Radon-Nikodym derivative: 

'^\Qt = £{- 0^'^-w^)^£:((^ - i)MQ)^i^f^, 

where M'^ := Nt — /q(1 — Nu)\udu is the compensated Q-martingale associated with A'", the first 
two Doleans-Dade exponentials are defined in (3.7) and (3.8), and 



n 

i=l 



exp 



Consequently, on top of the mark-to-market risk and event risk premia, the investor can potentially 
disagree with the market over jump intensity and jump size distribution of X', allowing for a richer 
structure of price discrepancy as well as the optimal liquidation strategy. 

As in Theorem 3.4, we compute the drift function in terms of pre-default price C and default 
risk premia, namely, 

G^it, x) = - (V,.C(t, x))^S(t, x)</.^'^(f, x) + {R{t, x) - C{t, x)) x) - x)) A(t, x) 

+ Ef / (C(t,x + ye,)-C(t,x))(A*(t,x)f(y)-A^(t,x)r(y))dyY (4.17) 
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where Cj := (0, . . . , 1, . . . , 0)^. We observe that the first two components of share the same 



functional form as G in (3.19), though the price function C is derived from the jump-diffusion 
model. Even if the investor and the market assign the same mark-to-market risk and event risk 
premia, discrepancy over jump intensity and distribution will yield different liquidation strategies. 
Under quite general affine jump- diffusion models, Duffie et al. [17J provide an analytical treatment 
of transform analysis, which can be used for the computation of our drift function. 

5 Optimal Liquidation of Credit Default Index Swaps 

We proceed to discuss the optimal liquidation of multi-name credit derivatives. In the literature, 
there exist many proposed models for modeling multiple defaults and pricing multi-name credit 
derivatives. Within the intensity-based framework, one popular approach is to model each default 
time by the first jump of a doubly-stochastic process. The dependence among defaults can be in- 
corporated via some common stochastic factors. This well-known bottom-up valuation framework 
has been studied in jl6| I38j. among many others. 

As a popular alternative, the top-down approach describes directly the dynamics of the cu- 
mulative credit portfolio loss, without detailed references to the constituent single names. Some 
examples of top-down models include O El [131 [361 EZ] ■ Iii particular, Errais et al. [20] proposed 
affine point processes for portfolio loss with self-exciting property to capture default clustering. 
For our analysis, rather than proposing a new multi-name credit risk model, we adopt the self- 
exciting top-down model from [20]. Also, we will focus on the optimal liquidation of a credit 
default index swap. 

First, we model successive default arrivals by a counting process {Nt)o<t<T, and the accu- 
mulated portfolio loss by T( = /i + . . . + /jvj, with each In representing the random loss at the 
nth default. Under the historical measure P, the default intensity evolves according to the jump- 
diffusion: 



d\t = k{e - \t)dt + a\J XtdW; + T] dTt, (5.1) 

where is a standard P-Brownian motion. We assume that the random losses (/„) are inde- 
pendent with an identical probability density function m on (0,oo). According to the last term 



in (5.1), each default arrival will increase default intensity A by the loss at default scaled by the 
positive parameter r/. This term captures default clustering observed in the multi-name credit 
derivatives as pointed out in ^20j. We assume a constant risk- free interest rate r for simplicity, 
and denote {'Ht)o<t<T to be the full filtration generated by N, T, and W^. 

The market measure Q is characterized by several key components. First, the market's mark- 
to-market risk premium is assumed to be of the form 

^Q,F ^ k{e - Xt) - K{e - Xt) 2) 

crVXt 



such that the default intensity in (5.1) preserves mean- reverting dynamics with different param- 
eters K and 9 under the market measure Q (see [H [27] for similar specifications). Secondly, we 
assume that the Q-default intensity is Xt := f^Xt, with a positive constant event risk premium. 
Thirdly, the distribution of random losses can be scaled under Q. Specifically, we assume that 
under Q the losses (In) admit the pdf m{z) := h{z)m{z), for some strictly positive function h 
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with h{z)'m{z)dz = 1 . Then, the Radon-Nikodym derivative associated with Q and P is 

'^^\nt = £{-ct>^^^w^)^kr, (5.3) 

where £{ - <j)^^^W^) is defined in and 

(rt poo \ 
J J {1- fih{z))Xum{z)dzdujYl{nh{li)). (5.4) 

Under the market pricing measure Q, the Q-default intensity evolves according to: 

dXt = K{fie - Xt)dt + ay^t dW^ + fir] dTt, (5.5) 

where := t^f + £ is a standard Q-Brownian motion. Similarly, we can define the 

investor's pricing measure Q through the investor's mark-to-market risk premium (p'^'^ as in (5.2) 
with parameters k and 9; default intensity Xt = jlXt with constant event risk premium jl; and loss 
scaling function h so that the loss pdf m{z) = h{z)rh{z). 

The credit default index swap is written on a standardized portfolio of H reference entities, 
such as single-name default swaps, with same notational normalized to 1 and same maturity T. 
The investor is a protection buyer who pays at the premium rate in return for default payments 
over (0, T]. Here, the default payment is assumed to be paid at the time when default occurs, 
and the premium payment is paid continuously with premium notational equal to H — Nt- 

The market's cumulative value of the credit default index swap for the protection buyer is 
equal to the difference between the market values of the default payment leg and premium leg, 
namely, 

pCDX ^^Q^ f e-"("~*) dT„|7^t} -E«3{p[J^ / e'^'^^'-'^H - Nu)du\'Ht}, t<T. (5.6) 

J{0,T] J{0,T] 



Hence, similar to (2.7), the protection buyer solves the following optimal stopping problem: 

y^CD^ = esssup E^{e-'-(--*)pf \nt}. (5.7) 

The associated delayed liquidation premium is defined by 

^CDX ^yCDX _pCDX_ (5 8) 

The derivation of the optimal liquidation strategy involves computing the market's ex-dividend 
value, defined by 

^CDX f g-r(«-i) _ jgQ r m f ^-r{u-t) ^jj _ |^ 1 ^ (5 9) 

J{t,T] J{t,T] 



Proposition 5.1. The market's ex-dividend value of the credit default index swap in (5.9) can 
be expressed as Cf^^ = C^^^{t,Xt,Nt), where 



C^^^{t, X, n) = k2{t, T)X + ki{t, T)n + A;o(t, T), (5.10) 
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for t <T, with coefficients 

Ut,T) = {cr+p^){ , , . + 

A:i(t,r) 



1 



p(/) + r) 
p^(l-e-^(^-*)) 



ce 



-r(T-i) 



/co(t,r) 



/ 

P V 



pr r(/3 + r)^ 

e-p{T-t) T-t 1 



(l_e-/'(^-*)), (5.11) 



(5.12) 



p{p + r) 



r/>' r2(r + p)- 



+ ce' 
and constants 



„-p(T-t) 1 



zm{z)dz, 



(1 



-r(T-t) 



and p = K — prjc. 



(5.13) 



(5.14) 







Proof. Using integration by parts, we re- write the market's ex-dividend value as 

(5.15) 

Hence, the computation of C*^^^ involves calculating E^^jA'^u | %t} and E^^jT^ | u > t. Since 
default intensity A follows a square-root jump-diffusion dynamics, these conditional expectation 
admit the closed-form formulas (see e.g. Section 4.3 of |20j): 



E^{A^„ \\t = \Nt = n, Tt = v] = A{t, u) + B{t, u)X + n, 
E'Q{T„ I At = A, iVt = n, = i;} = c^(t, u) + cB{t, u)X + v, 



(5.16) 
(5.17) 



for t < u < T, where 



A{t,u) 



Kp6 



-{K—firic)(u—t) 2 



K — pr]c K — prjc 

B{t,u) = (1 -e-^^-^''^)^"-*)) 

K — prjc 



+ u-t), 



(5.18) 
(5.19) 

Here, c is the market's expected loss at default given in (5.14). Substituting (5.16) and (5.17) 
into (5.15), we obtain the closed-form formula for market's ex-dividend value in (5.10). □ 

As a result, the ex-dividend value C'-'^^ is linear in the default intensity A^ and number of 
defaults Nf. Next, we characterize the optimal corresponding liquidation premium and strategy. 



Theorem 5.2. Under the top-down credit risk model in (5.1), the delayed liquidation premium 
associated with the credit default index swap is given by 



I e-'-(«-*)G^^^(n,A„)dn|At = A}, 



(5.20) 
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where 



A) = ( {f^7jk2{t,T) + l)(^ - c) + h{t,T){f^ - 1) - h{t,T){~K - ) A 
+ k2{t,T)fi{ke - kO), 



(5.21) 



with c := z'm{z)dz. If G'-''^^ {t, X) > V(t,A), then it is optimal to delay the liquidation till 
maturity T. If G'-''^^ [t, X) < V(t, A), then it is optimal to sell immediately. 



Proof. In view of the definition of L^^^ in (5.8), we consider the dynamics of p^"^^. First, it 
follows from (5.6) and (5.9) that 



^.riu-t)pCDX ^ ^-r{u-t)^CDX ^ f ^-r{v-t) ^^^^ -p^{H- K)dv) 



(5.22) 



Using (5.22) and the fact that e '"*pC'DX jg Q-martingale (whose SDE must have no drift), we 
apply Ito's lemma to get 



^-r{T-t)pCDX _ pCDX 



fifiCDX 



+ 



+ 



V e-^("-*)(T„-T„j- r r 

Jt Jo 



zm{z)Xu dzdu 



t<U<T 



e-^("-*) (C^«^(n, A„, 7V„) - C^^^(n, A„_, iV„_)) 



t<U<T 



t JO 



g-r(«-t) ^c^DX^^^ ^ ^^^^ + 1) - C^^^iu, A„, Nu))m{z)Xu dzdu 



(5.23) 



-r(u~t)(^ Ar„)^y;:A:(iT^„^ + G^^^(n, A„, iV„)dn) 



5A 



+ 



+ 



t<U<T 



-r{u—t) 



(T.-T„_) 



t JO 



zm(z)A„ dzdu 



^ (cC^^(^, A., iV„) _ c^^^(n, A._, iV._)) 

^-r{u-t) ^c^DX^^^ ^ ^^^^ + 1) - C^^^(u, A„, iV„))m(z)A„ dzdu 



t<U<T 
T fOO 



It JO 

for t < T < T, where 

f)f<CDX 

G^^^(t,X,n) := — (t, X,n)((kd - K9)fi - (k - k)X) 

oX ' 



(5.24) 



(5.25) 



+ / (z + C^'^^(t,A + /xr/z,n + l)-C^'^^(t,A,n))(^m(z)-m(z))Adz. 



Note that the two compensated Q-martingale terms in (5.23) account for, respectively, losses and 
changes in C^^^ value due to default arrivals. The second equation ( 5.24[ ) follows from change 
of measure from Q to Q . 

By Proposition 5j] the terms ^^^^^ and C*^^^ {t, X + ^r]z,n + 1) - C^^^ {t, X,n) do not 
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depend on n. Consequently, G^^-^ does not depend on n, and admits the closed-form formula 



(5.21) upon a substitutfon of (5.10) into (5.25) 



By taking the expectation on both sides of (5.24) under Q, the delayed liquidation premium 



l^CDX gai^ig£gg (5.20) and depends only on t and A. If G'^^^ > 0, then the integrand in (5.20) 



is positive a.s. and therefore the largest possible stopping time T is optimal. If G^^-^ < 0, then 
r* = t is optimal and Lf'^^ = a.s. □ 



We observe that the drift function consists of two components. The first component in (5.25) 



accounts for the disagreement between investor and market on the fluctuation of market ex- 
dividend value, while the second integral term reflects the disagreement on the jumps of market's 
cumulative value arising from the losses at default and the jumps in the ex-dividend value. Even 
though the market's cumulative value p^^^ in (5.6) and the optimal expected liquidation value 
V^^^ in Kt\ are path-dependent, both the delayed liquidation premium L^^^ in KM and 



To obtain the variational inequality of L'^^^ ^ we recall that A 
of default intensity A: 



qCDX jj-^ (5.21) depend only on t and A due to the special structure of C^^-^ given in (5.10) 

jlX/^i and the 



-dynamics 



dXt = k{fi6 — Xt)dt + (T \/ fiXt dW^ + ^rjdTt- 



The delayed liquidation premium as a function of time t and Q-default intensity A 

satisfies the variational inequality 



mm 



- k{fie - A) — ^ ;^T^ + rL^^^ 



dt 
/iA 
^ io 



dX 

(L^^^{t,X + fMr]z) 



2 aA2 

L^^^{t,X))m{z)dz 



G 



CDX 



CDX 



0, (5.26) 



for (t. A) G [0, T) X M, with terminal condition L{T, A) = for A G M. 

We consider a numerical example for an index swap with constant losses at default. In this 
case, the integral term in (5.26) reduces to L^^^ {t, X + jirjc) — L'^^^ {t, A), where c is the constant 



loss. We employ the standard implicit PSOR iterative algorithm to solve L^^^ by finite difference 
method with Neumann condition applied on the intensity boundary. There exist many alternative 
numerical methods to solve variational inequality with an integral term (see, among others, [2lll2j). 
We apply a second-order Taylor approximation to the difference L^^-^{t, X + firjc) — L'-'^^ [t, A) ~ 
dxL^^^ (t, X)fi7]c+ldx\L'^^^ {t, X){fn]c)^ . In turn, these new partial derivatives are incorporated 



in the existing partial derivatives in (5.26), rendering the variational inequality completely linear 



in A, and thus, allowing for rapid computation. 

We denote the investor's sell region S and delay region T> by 



S^^^ = {it,X) G [0,T] X 

v^^^ = {{t,x) G [o,r] X 



L^^^{t,X) = 0}, 
L^^^{t,X) > 0}. 



(5.27) 
(5.28) 



On the other hand, we observe from (5.10) a one-to-one correspondence between the market's 



ex-dividend value C'~'^^ of an index swap and its default intensity A for any fixed t < T, namely, 

C^'^^ - ki{t,T)n-ko{t,T) 



X 



k2{t,T) 



(5.29) 



Substituting (5.29) into (5.27) and (5.28), we can describe the sell region and delay region in 
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terms of the observable market ex-dividend value C'-^^^ . 

In Figure |4j we assume that the investor agrees with the market on all parameters except the 
speed of mean reversion for default intensity. In the case with k = 0.5 < 1 = k (Figure [4] (left)), 
the investor's optimal liquidation strategy is to sell as soon as the market ex-dividend value of 
index swap C^^^ reaches an upper boundary. In the case with k = 1 > 0.5 = R (Figure |4] 
(right)), the sell region is below the continuation region. 




012345 012345 
time time 



Figure 4: Optimal liquidation boundary in terms of market ex-dividend value of an index swap. We take 
T = 5, r = 0.03, = 10, = 0.25, a = 0.5, c = 5 = 0.5, pg" = 0.02, 6* = ^ = 1, and /i = fi = 1.1. 
Left panel: When k — 0.5 < 1 = k, liquidation occurs at an upper boundary that decreases from 3 to 
over t e [0,5]. Right panel: When k = 1 > 0.5 = k, the index swap is liquidated at a lower liquidation 
boundary, which decreases from 1.9 to over time. In both cases, the dashed line defined by G'-"^-^ — 
lies within the continuation region. 

In summary, we have analyzed the optimal liquidation of a credit default index swap under 
a top-down credit risk model. The selected model and contract specification give us tractable 
analytical results that are amenable for numerical computation. The top-down model implies 
that the underlying credit portfolio can experience countably many defaults. As argued by Errais 
et al. [20], this feature is innocuous for a large diversified portfolio in practice since the likelihood 
of total default is negligible. 

Our analysis here can be extended to the liquidation of CDOs. Consider a tranche with lower 
and higher attachment points Ki^K2 E [0, 1] of a CDO with H names and unit notionals. The 
tranche loss is a function of the accumulated loss T^, given by Lt = (T^ — KiH)^ — (Tt — K2H)^, 
t G [0,T]. With premium rate p™, the ex-dividend market price of the CDO tranche for the 
protection buyer is 

C^^O(t,At,Tt)=E^{ / e-'-^"'*) dLjl^ J - E«2{p- / e-''~''''\H{K2 - K,) - L^) du\nt} . 

J{t,T] J{t,T] 

Hence, the CDO price is a function of the accumulated loss T, as opposed to N in the case of 



CDX (see Proposition 5.1 above 



6 Optimal Buying and Selling 

Next, we adapt our model to study the optimal buying and selling problem. Consider an investor 
whose objective is to maximize the revenue through a buy/sell transaction of a defaultable claim 



{Y, A, R,Td) with market price process P in (2.4). The problem is studied separately under two 



scenarios, namely, when the short sale of the defaultable claim is permitted or prohibited. We 
shall analyze these problems under the Markovian credit risk model in Section [3] 
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If the investor seeks to purchase a defaultable claim from the market, the optimal purchase 
timing problem and the associated delayed purchase premium can be defined as: 

= essM E^{e--ft'' ''-'^"P^tlGt], and := - > 0. (6.1) 
6.1 Optimal Timing with Short Sale Possibility 

When short sale is permitted, there is no restriction on the ordering of purchase time r'' and sale 
time T'^. The investor's investment timing is found from the optimal double-stopping problem: 

Ut:= esssup E^je-^'t"' '^^'^^P^. - e-^i^'^^^^P^^igJ. 

Since the defaultable claim will mature at T, we interpret the choice of r** = T or = T as no 
buy/sell transaction at T. 

In fact, we can separate U into two optimal (single) stopping problems. Precisely, we have 

Ut = (esssup E^{e-^t'''-'^''Prs\gt} -Pt) + {Pt - essinf {e' '■^'^'^ P^tlGt}) 

= Lt + L\. (6.2) 

Hence, we have separated lA into a sum of the delayed liquidation premium and the delayed 
purchase premium. As a result, the optimal sale time r^* does not depend on the choice of the 
optimal purchase time r^*. 

The timing decision again depends crucially on the sub/super-martingale properties of dis- 
counted market price under measure Q. Under the Markovian credit risk model in Section 3, we 
can apply Theorem 3.4 to describe the optimal purchase and sale strategies in terms of the drift 



function G(t, x) in (3.19). 



Proposition 6.1. If G{t,x) > V(t,x) G [0, T] x M", then it is optimal to immediately buy 
the defaultable claim and hold it till maturity T, i.e. r^* = t and r** = T are optimal for Ut- 
If G{t,x.) < V(t,x) G [0,T] X M", then it is optimal to immediately short sell the claim and 
maintain the position till T, i.e. r** = t and t^* = T are optimal for Ut- 



6.2 Sequential Buying and Selling 

Prohibiting the short sale of defaultable claims implies the ordering: < r'^ < T. Therefore, the 
investor's value function is 

Ut:= esssup E"^ {e- Pr^ - e- ''^'^^ P^^lGt} . (6.3) 

The difference lit — Ut > Q can be viewed as the cost of the short sale constraint to the investor. 
As in Section^ we adopt the Markovian credit risk model, and derive from the Q-dynamics 
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of discounted market price in (3.20) to obtain 



Ut = ess sup 



L{t<rd} ess sup 



G{u,^u)du\Ft} ■ 



3.4) 



Using this probabilistic representation, we immediately deduce the optimal buy/sell strategy in 
the extreme cases analogues to Theorem |3.4[ 



Proposition 6.2. IfG{t, x) > \/{t, x) € [0, T] xM", then it is optimal to purchase the defaultahle 
claim immediately and hold until maturity, i.e. t^* = t and t^* = T are optimal for Ut- 

If G{t,x) < V(t,x) G [0,7"] X M", then it is optimal to never purchase the claim, i.e. 
^b* — ^s* _ rp optimal for Ut . 

Define U{t,'Kt) as the pre-default value of Ut, satisfying Ut := l{t<T-^}C/(t, Xj). We may view 
U{t,Xt) sequential optimal stopping problem. 



Proposition 6.3. The value function Ut in (6.3) can be expressed in terms of the delayed liqui- 
dation premium L in (3.21). Precisely, we have 



C/(t,Xt) = esssup E^{e--ft (""+^")'^%.| Ji}. 
T''ert,T 



(6.5) 



Proof. We note that, after any purchase time r^, the investor will face the liquidation problem 
V^b in (2.7). Then using repeated conditioning, Ut in (6.3) satisfies 

Ut = ess sup (e- ^t' ^-d^E^{e- T"' ^'"'^"P^. \g^b } - e~ ^""^"P^;,) l^t } (6.6) 

T^eTt,T,T^^T^b^T 

< esssup E^{e- ^-'^''{y^b - P^b)\Qt] (6.7) 

= esssup E^je-ii^'^-^^L^bigt} = l|t<,j esssup E^{e-i't'''(""+^")'^"4i,| J^}. (6.8) 

On the other hand, on the RHS of we see that V ^b = E'Qje" i'.'"'' '■"^"P^.. I^^b}, with the 

optimal stopping time r** := inf{t > r** : Vj = Pj} (see (2.10)). This is equivalent to taking the 
admissible stopping time r'^* for Ut in (6.3), so the reverse of inequality (6.7) also holds. Finally, 
equating (6.6) and (6.8) and removing the default indicator, we arrive at (6.5). □ 



According to Proposition 6.3, the investor, who anticipates to liquidate the defautable claim 
after purchase, seeks to maximize the delayed liquidation premium when deciding to buy the 
derivative from the market. The practical implication of representation (6.5) is that we first solve 
for the pre-default delayed liquidation premium L(f,x) by variational inequality (3.22). Then, 
using L{t,x) as input, we solve IJ by 



mm 



dU 



t, x) - Cr -M{t, x), U{t, x) - L{t, x) = 0, (t, x) e [0, T) x 



3.9) 
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where £^ ^ is defined in (3.17), and the terminal condition is [/(T, x) = 0, for x G M". In other 
words, the solution for L{t,x) provides the investor's optimal liquidation boundary after the 
purchase, and the variational inequality for U{t,x.) in (6.9) gives the investor's optimal purchase 
boundary. 

In Figure [5} we show a numerical example for a defaultable zero-coupon zero-recovery bond 
where interest rate is constant and A follows the CIR dynamics. The investor agrees with the 
market on all parameters except the speed of mean reversion for default intensity. When kx < k\, 
the optimal strategy is to buy as soon as the price enters the purchase region and subsequently sell 
at the (higher) optimal liquidation boundary. When k\ > R\, the optimal liquidation boundary 
is below the purchase boundary. However, it is possible that the investor buys at a lower price 
and subsequently sells at a higher price since both boundaries are increasing. It is also possible 
to buy-high-sell-low, realizing a loss on these sample paths. On average, the optimal sequential 
buying and selling strategy enables the investor to profit from the price discrepancy. Finally, 
when short sale is allowed, the investor's strategy follows the corresponding boundaries without 
the buy-first/sell-later constraint. 




Figure 5: Optimal purchase and liquidation boundaries in the CIR model. The common parameters are 
T = 1, r = 0.03, cr = 0.07, = /i = 2, and 6'a = = 0-015. Left panel: When ka = 0.2 < 0.3 = ka, 
the short sale constraint moves the purchase boundary higher. Both purchase boundaries, with or without 
short sale, are dominated by the liquidation boundary. Right panel: When ka = 0.3 > 0.2 = ka, the short 
sale constraint moves the purchase boundary lower. The liquidation boundary lies below both purchase 
boundaries. 

Remark 6.4. In a related study, Leung and Ludkovski [32j also discuss the problem of sequential 
buying and selling of equity options without short sale possibility and with constant interest 
rate. In particular, the underlying stock admits a local default intensity modeled by X{t,St), a 
deterministic function of time t and current stock price St- In contrast, our current model assumes 
stochastic default intensity At = A(t,Xt) and interest rate r(t,Xt), driven by a stochastic factor 
vector X. Hence, our optimal stopping value functions and buying/selling strategies depend on 
the stochastic factor X, rather than the stock alone as in |32j . 

7 Conclusions 

In summary, we have provided a flexible mathematical model for the optimal liquidation of various 
credit derivatives under price discrepancy. We have identified the situations where the optimal 
timing is trivial and also solved for the cases when sophisticated strategies are involved. The 
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optimal liquidation framework enables investors to quantify their views on default risk, extract 
profit from price discrepancy, and perform more effective risk management. Our model can also 
be modified and extended to incorporate single or multiple buying and selling decisions. 

For future research, a natural direction is to consider credit derivatives trading under other 
default risk models. For multi-name credit derivatives, in contrast to the top-down approach taken 
in Section 5, one can consider the optimal liquidation problem under the bottom- up framework. 
Liquidation problems are also important for derivatives portfolios in general. To this end, the 
structure of dependency between multiple risk factors is crucial in modeling price dynamics. 
Moreover, it is both practically and mathematically interesting to allow for partial or sequential 
liquidation (see e.g. [MIES])- On the other hand, market participants' pricing rules may vary due 
to different risk preferences. This leads to the interesting question of how risk aversion influences 
their derivatives purchase/liquidation timing (see e.g. [33] for the case of exponential utility). 
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